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GRID SPACING CONTROL WITH VARIATION DIMINISHING SPLINES 

By 

Philip W. Smith* 

The purpose of this research is to improve methods used 
to specify and control two and three dimensional grids on 
which numerical solutions of partial differential equations 
may be obtained. While initially focusing on grid generation, 
we anticipate that it will evolve into a consideration of 
the interaction of grid generation with the solution of a 
partial differential equation. 

The research during this period concentrated on using 
the multi-surface method of grid . aeration, developed by 
Peter Eiseman , to continuously patch a grid onto em existing 
grid. In the resulting grid the r laments of the Jacobian 
matrix must be continuous across the boundary between the 
original grid and the patched grid. 

This research is motivated by the desire of engineers 
to account for the affects of wind tunnel walls in physical 
experiments by comparing numerical solutions using an original 
grid whose boundary represents the wind tunnel walls with 
numerical solutions using a much larger grid formed by 
patching a grid around the original with exterior free stream 
conditions. 

Programs were written which accept as input the coordinates 
of the original grid and the desired new boundary, and then 
use the three-surface or four-surface version of the multi- 
surface method to extend the original grid out to a new 
boundary . 

*Professor, Department of Mathematical Sciences, Old Dominion University, 
Norfolk, Virginia 23508. 
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The multi-surface transformation generates coordinates 
between an inner boundary surface ^^^(s) and an outer boundary 
surface 5jj(s). Intermediate surfaces , . . . 

introduced to control the coordinate representation between 
?^(s) and Each surface representation is of the form 



X, (s)‘ 

Yj^(s) 


, k= 1,2, . . . ,N 


0 ^ Sj^ ^ 1 

i= 1,2 


The physical domain can be written 

’x(Pj^(s) , . . . »Ti)" 

P(s,n) = |y(P.j^(s) ,...,Pjj(s) ,n) 

z(5j^(s) ,««>,^^(s) ,r|) 


where ^ 1 for i= 1,2 and 0 s: n ^ 1. The variable n is 

the independent variable which spans between surfaces. 

The multi-surface method is explained in more detail in 
"Mesh Generation Using Algebraic Techniques" by Peter Eiseman 
and Robert Smith which can be found in Numerical Grid Generation 
Techniques , NASA Conference Publication 2166. 

The intermediate surfaces are not coordinate surfaces. 


but are surfaces used to determine a field of tangent vectors 
to the coordinate curve spanning across the surface. The 
vector field of tangents is given by 


Vj^(s) = Aj^(Pj^^^{s) - Pj^(s)) k= 1,2,...,N-1 


where the Aj^ are scalars determining only the magnitude of 


the vectors 
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Using the independent variable n .for the spanning direction 
we can choose ^ 2 ^ ••• < to correspond to the 

A sufficiently smooth vector field V(rifS) is obtained by 
a smooth interpolation V(t1j^,s) = Vj^ (s) . The first derivative 
of P{s,Ti) with respect to n is given by 


1 ) 




N-1 

^ where 

k=l ^ ^ 


= { 


1 

0 


if j = k 
if j?^k • 


Integrating 1) with an initial value of P 2 ^(s) yielas 

N-1 

p(s ,n) = P^(s) + 

where 


G^(n) = 


ii)j^(x)dx. 


If the interpolants are polynomials in n then the coordinate 
curves which cannect the bounding surfaces are defined by 
polynomials n of one degree greater. So in particular, using 
the three-surface method with 0, ti2= 1/ (n) = 1 “ n and 

4^2 tn) = n we get 

2 

2) Gj^(n) = n- ^ 


\ 
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Also, note that if the Aj^ are chosen so that 
is 1, then ^ (s,Tijg_jL) = . 

We then have 

p{s,n) = + V^>> 

which is Eiseman's three-surface transformation. 

Substituting 2) for Gj^, k= 1,2 and rewriting P(s,n) yields 

P(s,n) = I B. (n)P>{s) where 
k=l ^ ^ 

(n) = 1 - 2n + 

82 (n) = 2n - 2n^ 

BjCti) = . 

Using the four-surface method with 0, ^ 2 ~ ’^ 3 “ 

ti)^(n)= ( 1 - n) (j- n) f t |/2 (n) = n d - n) , and t|; 2 (ri)= (n-j)ri yields 

Gl(n) = 

G^cn) = 

03(11) = . 

So the four-surface transformation is given by 

3 G. ( ri ) ^ ^ ^ ^ 

P(S,T,) = ■ 
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After substitution and rewriting we have 

?(s,n) = I D, (ri)?. (s) where 
k»l ^ ^ 

Dj^(n) = 1- 6n+ 9n^- 4n^ 

D 2 (n) = 6n- 12n^+ 6n^ 

D3(n) = 6n^“ 6n^ 

D^Cn) = -3n^+ 4n^ 

In order to use the three-surface method to patch a grid 
around an existing grid we choose to be the next to the 

last row of the original grid and to be the boundary of 

the new grid. The vector s becomes a scalar variable t, 

0 s t s 1, since we are dealing with curves. To obtain 
go through the following steps: 

d(l.(t) 

Compute rift(t) = — r rr~^ where 

° d(?^(t),?3(t)) 

d(A,$) = distance from A to 1. 

Denote the last row of the original grid by 

Pick P, (t) to solve 
^ 3 

S(t) = z B. (n«(t))P. (t) . 
k=l ^ ^ ^ 

We will then have P (t,0) = Pj^(t) , ^{t, Tig (t) )* ^(t.) , and 
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Along each connecting curve we choose m equally 
spaced values / i* 0,...,ni~l from nQ(t) to 1 with ra= 1/t1q(0) 
and compute the corresponding physical points P(t,n£)* ^yJt^Ti^’)^ 
We then make a table consisting of accumulated arclengths 
alpha,, i= 0,...,m-l along C. from the point where n = nn(t) 

X ^ w 

to where n = 1 versus the corresponding value Q"f ti. 

Using the arclength, ARC, from Pj^(t) to Q(t) we use the 
Newton-Raphson method to solve 

arc = l)/(e^^^^ - 1) 

for k(t) where An = l/(m-l). 

Next we substitute m values of n from 0 to 1 in increments 

of An into l)/(e^^^^ - 1) to obtain an exponential 

distribution of arclength along C,^ from Pj^(t) to P^Ct). This 

is achieved by interpolating the values from the exponential 

expression into the table created above to find the appropriate 

n*(t,n) corresponding to a particular arclength. The value 

of n*(t,n) is then substituted into P(t,n) to obtain the 

physical coordinates for the patched grid. Our transformation 

can then be written 

^ ^ 3 

F(t,n) = P(t,n*) = I B, Cn* (t,n) ]Pv (t) 

k=i ^ ^ 

where 0^nsl,0:st<l. 

The function n*(t,n) is defined by n* (t,n) = u (a) where 
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u(a) 


alpha^^j^) + 
ni_2^(alphai~ a) ]/ 
(alpha^- alpha^_^) 


for a = alpha^ 


for alpha^_j^< a< alpha^. 


This transformation guarantees continuity in the spacing 
between rows as one moves from the original grid to the 
patched grid, but continuity in the spacing between columns 
appears to be dependent on the distribution of points along 
the new boundary. Our current research centers on finding 
a reliable method for choosing the best boundary points. 
Hopefully, an interactive method for choosing the points 
can be inserted into the current progrcuns. 

All progrcuns were done at NASA Langley Research Center 
on the Interactive Processing System using a Tektronix 4014 
graphics terminal and a Prime 750 computer. This system 
allows one to interactively generate grids fairly easily. 



